Nanogranular Thermoelectrics 
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We investigate thermopower and thermoelectric coefficient of nano-granular materials at large 
tunneling conductance between the grains, 2> 1. We show that at intermediate temperatures, 
T > grS, where S is the mean energy level spacing for a single grain, electron-electron interaction 
leads to an increase of the thermopower with decreasing grain size. We discuss our results in the light 
of new types of thermoelectric materials and present the behavior of the figure of merit depending 
on system parameters. 
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The search for more efficient thermoelectric materials 
has had little success during the last several decades since 
bulk materials are limited in their performance by the 
Wiedemann-Franz law that connects electric to thermal 
conductivity in such a way as to defeat all attempts at im- 
proving the dimensionless figure of merit ZT = S'^aT / k, 
where S is the thermopower or Seebeck coefficient, tj is 
the electric conductivity, and k is the thermal conductiv- 
ity [H, [3, H, 0, [1] . To be competitive compared with con- 
ventional refrigerators, one must develop thermoelectric 
materials with ZT > 2. Although it is possible in prin- 
ciple to develop homogeneous materials with ZT > 2, 
there are no candidate materials on the horizon. Thus, 
one needs to search for inhomogeneous/granular thermo- 
electric materials in which one can directly control the 
system parameters. 

Most theoretical progress was archived by numerical 
solution of phenomenological models d, 0|- However, 
no analytical results obtained from a microscopic model 
for coupled nanodot/grain systems is available till now. 
Thus, the fundamental question that remains open is 
how thermoelectric coefficient and thermopower behave 
in nanogranular thermoelectric materials. Here, we make 
the first step towards answering this question for granu- 
lar metals at intermediate temperatures by generalizing 
our approach 0] recently developed for the description 
of electric @ and heat transport [l^l- In particular, we 
will answer the question to what extend quantum and 
confinement effects in nanostructures are important in 
changing ZT. 

In this paper we investigate the thermopower S, ther- 
moelectric coefficient rj, and the figure of merit ZT of 
granular samples focusing on the case of large tunnel- 
ing conductance between the grains, 3> 1. With- 
out Coulomb interaction the granular system would be a 
good metal in this limit and our task is to include charg- 
ing effects in the theory. 

Main results. — The main results of our work are 
as follows: (i) We derive the expression for the ther- 
moelectric coefficient 77 of granular metals that includes 
corrections due to Coulomb interaction at temperatures 



T > gxS, where S is the mean level spacing of a single 
grain 
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Here 77(0) = -{t:^ /3)egTa^-'^{T/eF) is the thermoelec- 
tric coefficient of granular materials in the absence of 
electron-electron interaction with e being the electron 
charge, a the size of a single grain, d = 2, 3 the dimen- 
sionality of a sample, ep being the Fermi energy, and 
Ec = e^/a is the charging energy. 

(ii) We obtain the expression for thermopower S of 
granular metals 



^ = 5(°)fl-^ln-9^^'= 
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where 5''^°-' = —{n^ /6){l/e){T /ep) is the thermopower of 
granular metals in the absence of Coulomb interaction. 
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FIG. 1: Plots of the dimensionless figure of merit Z/Z^ ' vs. 
grain size a (in nm)- for different values of dimensionless tun- 
neling conductance gr (see legend): the upper panel is for 3D 
and the lower for 2D. All curves are plotted for T = lOOK. At 
this temperature, the dimensionless bare figure of merit for 
granular metals is Z^°^T « lO""*. 
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(iii) Finally, wc find the figure of merit to be: 



over all grains in the system and 
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where Z^^'^T = {w'^ /l2){T/eF)'^ is the bare figure of merit 
of granular materials and 7 0.355 is a numerical coef- 
ficient. In Fig. [1] we plot Z/Z(°) as a function of the 
grain size a for different tunneling conductances gr at 
fixed temperature lOOK. We find, that the influence of 
granularity is most effective for small grain sizes and the 
presence of Coulomb interaction decreases the figure of 
merit. 

At this point we remark that all results are obtained 
in the absence of phonons which become relevant only at 
higher temperatures. At the end of this paper we will 
briefly discuss their influence. 

Our main results, Eqs. H]) - ([3]), are valid at interme- 
diate temperatures, T > gx^- At these temperatures the 
electronic motion is coherent within the grains, but co- 
herence does not extend to scales larger than the size a 
of a single grain Q. Under these conditions, the electric 
conductivity a and the electric thermal conductivity k 
are given by the expressions 0, 11 1 
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where ct(°) = 2e^gTa^-'^ and k'") = L^a^^'^T are the elec- 
tric (including spin) and thermal conductivities of gran- 
ular metals in the absence of Coulomb interaction with 
Lq = 7r^/3e^ being the Lorentz number. We mention 
that at temperature T > g^i the correction to the ther- 
moelectric coefficient, Eq. ([T]), has a TlnT dependence 
in both d = 2, 3 dimensions which is similar to the re- 
sult for the electric conductivity, Eq. (j4a| . having a InT 
dependence in all dimensions as well. 

Model. — Next we introduce our model and de- 
scribe the derivation of Eqs. (P)-®: We consider a 
d— dimensional array of metallic grains with Coulomb in- 
teraction between electrons. The motion of electrons in- 
side the grains is diffusive and they can tunnel from grain 
to grain. We assume that the sample would be a good 
metal in the absence of Coulomb interaction. However, 
we also assume that gx is still smaller than the grain 
conductance go, meaning that the granular structure is 
pronounced and the resistivity is controlled by tunneling 
between grains. 

Each grain is characterized by two energy scales: (i) 
the mean energy level spacing 5, and (ii) the charging 
energy Ec = e} j a (for a typical grain size of a w lOnm 
Ec is of the order of 2000K) and we assume that the 
condition (5 <C i?c is fulfilled. 

The system of coupled metallic grains is described by 
the Hamiltonian H = '^.^Hi, where the sum is taken 



k j^i 
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The first term in the right hand side (r. h. s.) of Eq. ([5]) 
describes the i-th isolated disordered grain, a|j,(ai,fc) are 
the creation (annihilation) operators for an electron in 
the state k and £,k = k'^ /2m — ji with ji being the chemi- 
cal potential. The second term describes the charging en- 
ergy, Cij is the capacitance matrix and fii = d\f,di^k 
is the number operator for electrons in the z-th grain. 
The last term is the tunnel Hamiltonian where tij are 
the tunnel matrix elements between grains i and j. 

Derivation of kinetic coefficients. — The kinetic coef- 
ficients: electric conductivity cr, the thermoelectric co- 
efficient 77, and the thermal conductivity k are related 
to the Matsubara response functions L^"'^' with a, /3 G 
{e,/i} ii,E2 



j(e) 



-(L(^'=V(e'T)) V{eV) - (L^^''^ /(eT^)) VT, 
-(£('='* V(er)) V{eV) - {L^''''^ /T^) VT. (6) 



Here j^*^^ (j *•'*•') is the electric (thermal) current 



and V the electrostatic potential. From Eq 
one finds that a = L^'"''>/T, rj = L^''^^/T^, 
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77 = L'-'^'-'/i", and 
S = -AF/AT = L^^'^yiTL^""')), where the re- 
sponse functions are given by Kubo formulas L^"^") — 

dre'"'"^(rj(")(T)j(/')(0)) 

with Tr being the time ordering operator for the currents 
with respect to the imaginary time r. Thus, to calculate 
the thermoelectric coefficient rj and thermopower S of 
the granular metals one has to know the explicit form of 
the electric j*^'^^ and thermal j^'*) currents. 

The electric current jj*^' through grain i is defined as 
i't^ — ^jJif — edhi/dt — ie[ni,H]. Straightforward 



calculations lead to j 



k,q 



ai,k). 









jk 




^ Jh 










a) 


b) 



FIG. 2: 



Vertices corresponding to the thermal current op- 



erator, Eqs. ([7]): vertex (a) corresponds to j^"'"' and (b) to 
rpj^g solid lines denote the propagator of electrons, the 
thick wavy line describes Coulomb interaction, the tunneling 
vertices are described by the circles, ci;„ = nT{2n + 1) and 
Qm = 2nmT are Fermionic and Bosonic Matsubara frequen- 
cies respectively (n, m e Z). 



3 



For granular metals the thermal current operator 
Ji'*'' = Sj Jij'' can be obtained as follows. The energy 
content of each grain changes as a function of time, such 
that dHi/dt = i[Hi,H]. Energy conservation requires 
that this energy flow to the other grains in the system, 
^''^ Calculating the commutator [Hi,H], 



, , . .(/t) 4/1,0) 
we obtam j\- = 7, 



' where 
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where {A; B}+ stands for the anti-commutator. The 
contribution is the heat current in the absence of 

electron-electron interaction, while the second term 
appears due to Coulomb interaction. Equations ([7]) im- 
ply that the thermal current operator must be associ- 
ated with two different vertices in diagram representa- 
tion. Fig. m 

For large tunneling conductance, the Matsubara ther- 
mal current - electric current correlator can be analyzed 
perturbatively in l/gr, using the diagrammatic tech- 
nique discussed in Ref. that we briefly outline be- 
low. The self-energy of the averaged single electron 
Green function has two contributions: The first contri- 
bution corresponds to scattering by impurities inside a 
single grain while the second is due to processes of scat- 
tering between the grains. The former results only in 
small renormalization of the mean free time which de- 
pends in general on the electron energy u as = 
Tq^ [1 + {d/2 — l)u;/eF] which is a result of the renor- 
malization of the density of states at the Fermi surface. 

The diffusion motion inside a single grain is given by 
the diffusion propagator Pq"^ = r^^iriil, where fii is the 
bosonic Matsubara frequency for the (internal) Coulomb 
interaction. The coordinate dependence in Vq is ne- 
glected since in the regime under consideration all char- 
acteristic energies are smaller than the Thouless energy. 
The complete diffusion propagator is given by ladder dia- 
grams resulting in the following expression 15^^(51^, q) = 
T'wd^^il + ^qS), where eq = 2,gTX]a(l ~ cosqa) with a 
being the lattice vectors. The same ladder diagrams 
describe the renormalized interaction vertex. The in- 
teraction vertex is used to obtain the polarization op- 
erator, that defines the effective dynamically screened 
Coulomb interaction for granular metals, V{fli,q) = 
2 [i?~^(q) + 4eq/ + eqS)] , where the charging en- 
ergy in d = 2, 3 is given by ^^(q) = 2^ = T!^- 

However, there is an important difference between cal- 
culations of thermoelectric coefficient 77 and thermopower 
S and the calculations of the electric a and thermal k con- 
ductivities in Eqs. (|4al I4bp . Indeed, to calculate a and k 




FIG. 3: Diagrams describing the thermoelectric coefficient of 
granular metals at temperatures T > qtS: diagram (a) cor- 
responds to rjo in Eq. Diagrams (b)-(d) describe first 
order corrections to the thermoelectric coefficient of granular 
metals due to electron-electron interaction. The solid lines 
denote the propagator of electrons, the wavy lines describe 
effective screened electron-electron propagator, and the (red) 
triangles describe the elastic interaction of electrons with im- 
purities. The tunneling vertices are described by the circles. 
The sum of the diagrams (b)-(d) results in the thermoelectric 
coefficient correction 77'^' given in Eq. (|8]). 



it was sufficient to approximate the tunneling matrix el- 
ement tpq as a constant t which is evaluated at the Fermi 
surface and neglect variations of tpq with energy which 
occur on the scale T/ep- However, this approximation is 
insufficient for calculations of thermoelectric coefficient 
?7 and thermopower S since the dominant contribution 
to these quantities vanishes due to particle-hole symme- 
try such that that both quantities are proportional to 
the small parameter T/ep. Since it is necessary to take 
into account terms of order of T/ep in order to obtain a 
nonzero result for 77 and S the corresponding expansions 
must be carried out to this order for all quantities which 
depend on energy: the density of states, the relaxation 
time, and the tunneling matrix element. For the latter we 

obtain the expression: ^^(6,6) = ''l + ll? 





In the absence of the electron-electron interaction the 
thermoelectric coefficient is represented by diagram (a) 
in Fig. 1. Straightforward calculations of this diagram 
lead to the result for 77^°^ given below Eq. 

First order interaction corrections to the thermoelec- 
tric coefficient are only generated by diagrams (b) and 
(c) in Fig. [21 resulting after summation over Fermionic 
and Bosonic frequencies and analytical continuation in 



^(1) 



271 gT 









T 



(8) 



where the q-integration goes over the d-dimensional 
sphere with radius 7r/a. Diagram (d) in Fig. Ogives only 
contributions to the thermoelectric coefficient of order 
(T/ep)'^ and higher. Integrating over q in Eq. ([S]) we 
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obtain the following expressions: 
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which lead to Eq. 

Discussions. — In the presence of interaction effects 
and not very low temperatures T > gr^i granular met- 
als behave differently from homogeneous disorder met- 
als. However, in the absence of interactions the result 
for 77^°^ below Eq. ^ coincides with the thermoelectric 

coefficient of homogeneous disordered metals, 77^o^„j = 
— (2/9)epi?(ror), with pp being the Fermi momentum. 
One can expect that at low temperatures, T < gxS, even 
in the presence of Coulomb interaction the behavior of 
thermoelectric coefficient and thermopower of granular 
metals is similar to the behavior of j]hom and Shom, how- 
ever this temperature range is beyond the scope of the 
present paper . Our results for thermopower ([2]) and fig- 
ure of merit ([3]) show that the influence of Coulomb in- 
teraction is most effective for small grains. decreases 
with the grain size which is a result of the delicate com- 
petition of the corrections of thermoelectric coefficient 
([T|) and the electric conductivity (gaj). In particular, if 
the numerical prefactor of the correction to ry would be 
slightly smaller, the sign of the correction to S would 
change. 

Above we only considered the electron contribution 
to figure of merit. At higher temperatures T > T*, 

where T* ^ ^ 9TCpf^/lpha is a characteristic tempera- 
ture with Ipfi and Cpf^ being the phonon scattering length 
and phonon velocity respectively [lo| . phonons will pro- 
vide an independent, additional contribution to thermal 
transport, Kph = T^lph/ c^h- However, the phonon contri- 
bution can be neglected for temperatures gT^ < T < T* . 
A detailed study of the infiuence of phonons at high tem- 
peratures, including room temperature, will be subject of 
a forthcoming work. 

So far, we ignored the fact that electron-electron in- 
teractions also renormalize the chemical potential fi. In 
general, this renormalization may affect the kinetic coef- 
ficients: the thermal current vertex. Fig [21 as well as the 
electron Green functions depend on /z. In particular one 
needs to replace V{eV) — > V{eV + fj,) in Eqs. ©. To first 
order in the interactions, the renormalization of fi only 
leads to corrections to diagram (a) in Fig. [3l As it can 
be easily shown, for this diagram the renormalization of 
the heat and electric current vertices is exactly canceled 
by the renormalization of the two electron propagators. 
Therefore, the renormalization of the chemical potential 
does not affect our results in the leading order. 

Finally, we remark that the bare figure of merit Z^'^^T 
for granular metals at > 1 and lOOK is of the order 
of only 10^**. Therefore these materials are not suit- 
able for solid-state refrigerators, but should be replaced 
by granular semiconductors with (7^ < 1. Therefore 



we conclude this paragraph discussing the dimension- 
less figure of merit ZT of granular materials at weak 
coupling between the grains, ^ 1- In this regime 
the electronic contribution to thermal conductivity 
of granular metals was recently investigated in Ref. Il3l . 
where it was shown that ~ g\T^ jE"^. In this regime 
the electric conductivity of gra nular metals obeys the 
law a ~ gT exp {-Ec/T) [i,IlH. However, an expres- 
sion for the thermoelectric coefficient in this region is 
not available yet, but recently it has been proposed, 
based on experiment, that nanostructured thermoelectric 
materials in the low coupling region (A£PbmSbTe2+m, 
BisTca/SbaTea, or CoSbg) [ij, E, [3, [B, El can have 
higher figures of merit than their bulk counterparts. 

In conclusion, we have investigated the thermoelectric 
coefficient and thermopower of granular nanomaterials 
in the limit of large tunneling conductance between the 
grains and temperatures T > gr^. We have shown to 
what extend quantum and confinement effects in granu- 
lar metals are important in changing ZT depending on 
system parameters. 
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